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We consider the effect of the RKKY interaction between magnetic impurities on the eiectron 
relaxation rates in a normal metal. The interplay between the RKKY interaction and the Kondo 
effect may result in a non-monotonic temperature dependence of the electron momentum relaxation 
rate, which determines the Drude conductivity. The electron phase relaxation rate, which determines 
the magnitude of the weak localization correction to the resistivity, is also a non-monotonic function 
of temperature. For this function, we find the dependence of the position of its maximum on the 
concentration of magnetic impurities. We also relate the electron energy relaxation rate to the 
excitation spectrum of the system of magnetic impurities. The energy relaxation determines the 
distribution function for the out-of-equilibrium electrons. Measurement of the electron distribution 
function thus may provide information about the excitations in the spin glass phase. 
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I. INTRODUCTION 

Electron transport in normal metals is known to be 
very sensitive to the presence of magnetic impurities in a 
metal. Scattering of conduction electrons off such impu- 
rities scrambles the electron spin. A tiny concentration 
of magnetic impurities results in an observable effect - 
low-temperature saturation of the phase relaxation ratei. 
Magnetic impurities apparently also facilitate the energy 
transfer between electrons 2 * 3 . At higher concentrations 1 , 
a minimum in the temperature dependence of the resis- 
tivity becomes evident, which is a manifestation of the 
Kondo effect. These three observations fit very well with 
a picture of uncorrelated magnetic impurities. 

Investigation of the resistivity at even higher concen- 
tration of magnetic impurities, reveals deviations from 
the picture of uncorrelated localized magnetic moments. 
The temperature dependence of the resistivity, in addi- 
tion to the aforementioned minimum, develops a maxi- 
mum 5 J 6 J 7 J 8j9 i ower temperatures. The electron phase 
relaxation rate was recently measuredi2*ii on AuFe al- 
loys with magnetic impurity (Fe) concentration ranging 
from 7.1 to 60 ppm. At the upper end of this range, 
correlations between the localized spins may become im- 
portant, as evidenced by the temperature dependence of 
the resistivitj«£*ii. 

In this paper we investigate the effect of correlation 
between spins of magnetic impurities on the electronic 
transport properties of a metal. Specifically, we study the 
Kondo contribution to the Drude resistivity, the weak- 
localization correction to the conductivity, and the elec- 
tron energy relaxation rate. 

The leading mechanism of correlations is known to 
be the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action between impurities. The interaction lifts the de- 
generacy in the excitation spectrum of impurities. In this 
respect, the effect of RKKY interaction is somewhat simi- 
lar to the effect of an external magnetic field which causes 
Zeeman splitting of the spin states. It is known that the 



magnetic field reduces the electron relaxation rates at low 
energies&iSii&i^, and also suppresses the Kondo effect. 
Unlike the uniform Zeeman splitting, however, interac- 
tion between the spins results in a broad spectrum of 
energies of collective spin states. Therefore, the quan- 
titative manifestations of the RKKY interaction in the 
electron transport are different from that of the Zeeman 
energy. 

The difficulty of the low-temperature electron trans- 
port problem is associated with the complexity of the 
spin glass state. In this paper, we study in detail the 
transport at relatively high temperatures (or electron en- 
ergy transfers, in the case of energy relaxation). We per- 
form analytical calculation using the method of the virial 
expansionist in the RKKY interaction between mag- 
netic impurities. The method is based on the following 
concept. Since the RKKY interaction decreases fast as 
a function of the distance between impurities, the impu- 
rities have to be close to each other for the interaction 
between them to compete with thermal smearing and to 
affect transport properties of conduction electrons. We 
perform the virial expansion to the second order in the 
density of magnetic impurities, which corresponds to ac- 
counting for the interaction within impurity pairs. 

Electron scattering off magnetic impurities contributes 
to the temperature dependence of the resistivity: due to 
the Kondo effect, the resistivity increases as the tem- 
perature is lowered. Similar to the Zeeman splitting, 
the RKKY interaction between magnetic impurities may 
stop the development of the Kondo effect. The inter- 
play between the RKKY interaction and the Kondo ef- 
fect leads to a maximum in the temperature dependence 
of the resistivity^^ If the characteristic temperature of 
the spin glass formation exceeds significantly the Kondo 
temperature (high concentration of magnetic impurities 
n s ), then the resistivity has a maximum at a tempera- 
ture in the region of applicability of the virial expansion. 
On the contrary, at small n s the maximum of resistivity 
occurs at zero temperature. 
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Magnetic impurities also affect the magneto-resistance 
in weak magnetic fields. This magneto-resistance is due 
to the weak localization (WL) effect42*i^ Being an in- 
terference phenomenon, the WL is limited by the elec- 
tron phase relaxation. We calculate the phase relaxation 
rate taking into account the RKKY interaction between 
magnetic impurities. At high impurity concentration, the 
phase relaxation rate has a maximum, similar to the max- 
imum in the resistivity. At low concentration, unlike the 
Drude resistivity, the phase relaxation retains a maxi- 
mum at a finite temperature (of the order of the Kondo 
temperature) . 

The RKKY interaction lifts the degeneracy of the mag- 
netic impurity states. Therefore an out-of-equilibrium 
electron may loose its energy by exciting the impurity 
spin degrees of freedom. The corresponding relaxation 
rate is a function of the transferred energy. This rate pro- 
vides information about the spin excitations spectrum. 
We derive the corresponding kinetic equation for the elec- 
tron distribution function. We also make specific predic- 
tions for the relaxation rate at sufficiently large energy 
transfers, which may be accounted for by the virial ex- 
pansion. 

Before we proceed, we emphasize that in this paper 
we consider the effect of the interaction between Kondo 
impurities only on kinetic properties (conductance, en- 
ergy relaxation) of electrons in a metal. The interac- 
tion between magnetic impurities may also affect ther- 
modynamic properties (such as heat capacity, suscepti- 
bility, superconducting transition temperature) of con- 
duction electrons >iSiiL22i2i An analysis of the thermody- 
namic properties of electrons in dilute Kondo alloys was 
performed earlier, see e.g. ref. l2fl 

The paper is organized as follows. In the next sec- 
tion wc introduce the model and discuss the effect of 
interaction between magnetic impurities on a spin cor- 
relation function. In Sec. Ill we perform calculations 
of the resistivity correction due to the electron scatter- 
ing off magnetic impurities. Section IV contains analysis 
of the WL correction to the conductivity. In Sec. V we 
derive the kinetic equation for the system of magnetic im- 
purities and conduction electrons and discuss the energy 
exchange rate between these two subsystems. Section VI 
contains discussions and conclusions. 



II. MODEL 

The scattering of conduction electrons off a magnetic 
impurity is described by the following Hamiltonian 

H c = JSa, (1) 

where S is the spin operator of a magnetic impurity and 
<t is the spin operator of a conduction electron repre- 
sented in terms of the Pauli matrices {<J X , a y , a z }. The 
exchange constant J is renormalized due to the Kondo 
effect and varies as a logarithmic function of energy e 
of conduction electrons. At temperature T higher than 



the Kondo temperature Ik, the exchange constant for 
thermal electrons, e ~ T, is given by 

J = -ln~ x ^, T»T K , (2) 

where v is the Fermi density of states per spin degree of 
freedom. 

In order to evaluate the effects of electron scattering 
off magnetic impurities on electron transport, we intro- 
duce the statistical averages of impurity spin compo- 
nents, (S a Sp}, with respect to thermodynamic states of 
the magnetic impurity system. We show that the electron 
transport properties in a metal with strong spin-orbit 
coupling are determined by the following spin correlator 

K(t) = (S(0)S(t)). (3) 

In general, the correlation function K (t) can be rewritten 
in terms of exact quantum states |£) of the system of 
magnetic impurities: 

K(w) = 2tt | (Z\S\Of 6(Ei - E e - w). (4) 

Here is the energy of state and p£ is the density 
matrix p^ oc exp(—E^/T). 

For a free magnetic impurity with spin S (all spin 
states are degenerate) the spin correlation function K (t) 
does not depend on time and its Fourier transform has 
the form: 

K 1 (u) = 2ttS(S + 1)S(uj), (5) 

Using the Fermi golden rule and Eq. (JBJ, we obtain the 
following expression for the electron scattering rate off 
magnetic impurities 

— = 2Tnsn s J 2 S{S + 1). (6) 

Here n s is the magnetic impurity concentration per vol- 
ume. The quantity r s is the mean free time of scattering 
off magnetic impurities. 

In metals, the leading interaction between magnetic 
moments of impurities is described by the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) mechanism. The corre- 
sponding Hamiltonian has the form: 

H RKKY = Y / V(n j )S i S J . (7) 

ij 

The magnitude of the RKKY interaction is given by the 
following expression 

V(r) = ^-co Si p, (8) 

where <p changes fast on the length scale of the Fermi 
wavelength Af ■ The interaction constant Vq may be rep- 
resented in terms of the exchange constant and electron 
density of states v. 

W = ^f^, Jo(r) = - ] , (9) 
2ir v m[vF/(riK)\ 
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The typical value of the RKKY interaction between two 
impurities separated by distance l/nl^ 3 (average distance 
between impurities) is 

2n s 



n s V (r 



-1/3 



■ku ln 2 (w F n s 1/3 /T K ) 



(10) 



Due to the randomness of the RKKY Hamiltonian 
-ffrtKKY, finding of states |£) and corresponding energies 
is a hardly possible task. At sufficiently low tem- 
perature a phase transition into a spin glass state may 
occur with extremely complicated structure of the wave 
functions |£) and energy spectrum^ A transition tem- 
perature is comparable with the typical energy of the 
interaction between magnetic impurities T sg . 

We focus our attention on the high temperature limit 
T 3> T sg and calculate the spin-spin correlation function 
K(t), Eq. ip]l. using the virial expansion method. In this 
method the interaction between magnetic impurities is 
taken into account only if the splitting of the spin states 
due to the interaction exceeds the system temperature; 
otherwise the interaction does not significantly change 
the spin correlation function. Therefore, the interaction 
is important only for magnetic impurities in clusters of 
the size ~ y/Vo/T. For a uniform distribution of mag- 
netic impurities in the metal, the probability for a forma- 
tion of such a cluster of k impurities scales as (T ag /T) k ~ 1 . 
We consider only clusters containing two (fc = 2) mag- 
netic impurities. 

The energy states |£) of two interacting spins are clas- 
sified by the total spin J (J — 0,1, ... , 2S) and its pro- 
jection M on a fixed direction: |£) = | J, M). The energy 
spectrum is given by 



J( J + 1) - 2S(S • 



il = V(r)ej (11) 



Ej = V(r 



and is degenerate with respect to the projection M. The 
spacing between levels with different J is proportional to 
the magnitude of the RKKY interaction, Eq. |SJ 

According to Eq. the corresponding spin correla- 
tion function of a magnetic impurity within distance r 
from another magnetic impurity may be represented in 
the form 



K 2 (t,V) = > ; Ajj, 



2S 

= £ 

J,J'=0 



2J+1 

'W) > 



y{i(ej-ej,)t-ej/T} 



where the statistical sum Z(V) is 



2.V 



Z(V) = ^(2J + l)e- yej / T , 



(12) 



(13a) 



j=o 



and the matrix elements Ajj> are 

Ajj. = ^2(JM\S\J'M')(J'M'\S\JM). (13b) 



The analytical form of the matrix elements is presented 
in Appendix A. We emphasize that Ajji are independent 
of the pair spin projection M . 



As we will see in the following sections, the electron 
transport properties are determined by the spin correla- 
tion function K(u>), averaged over configurations of mag- 
netic impurities. We calculate K(u>) within the virial 
approximation: 

K(u) = + J p(r) [K 2 (uj, V(r)) - K^cj)} d 3 r, 

(14) 

where p(r) is the probability density for two magnetic 
impurities to be at distance r; for uniform impurity dis- 
tribution p(r) = n s . Averaging over the relative posi- 
tion of two magnetic impurities in Eq. I|14l) can be per- 
formed in two steps. First, we make the substitution 
r 3 = Vo cos (p/(Ty), and then we perform integration over 
the fast varying phase ip. As the result, the spin correla- 
tor has the formii 



K(u) = Ki(u) +SK 2 (u>), 



(15) 



where 



^ U) = to T v S(S + l) J* [p ^ y) _ m] (16a) 



3 T 

2S 



P(uj,y) = J2 AjJ ' 



2J+1 e- yej 
S(S + l)Z(Ty) 



JJ'=0 

x5{uj-Ty[ej-ej,}). 



(16b) 



We notice that the spin correlation function Eq. i|15|) 
increases as the frequency uj decreases: 



8vr T sg ^ \ej - ej, 



K(co) 



exp 



UJEj 



(ej-ej,)T 



Z{u/{ej-ej,)) 



, w?0. (17) 



One may expect that due to the 1/w 2 behavior of the spin 
correlation function, Eq. I|17|l. the virial approximation 
breaks down even at T > T sg . Nevertheless, due to the 
property P(uj,0) = S(ui), see Eq. (|A4() . the integrand 
in Eq. (|16a|l has no singularity at y — 0, and the slow 
modes (oj < T) of the spin system do not affect electron 
transport. 

Equation (|T5)l supplemented with Eqs. ljl3"|) and {Tfj)| 

determines the spin correlation function at high temper- 
ature T ^> T sg . Below we use these equations to describe 
the effect of interaction between magnetic impurities on 
electron transport in metals. 



III. RESISTIVITY 

The conductivity of a metal with isotropic impurities 
may be calculated according to the standard rules of the 
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diagrammatic technique, 
corrections we have 



Disregarding the interference 



G R (e,p)G A (e,p) dedp 
Tcosh 2 e/2T (2tt) 4 ' 



Here 



G R > A {e,p) 



e - £(P) - ^ A (e) 



(18) 



(19) 



is the retarded or advanced Green function averaged over 
disorder, £(p) = «f(|p| — Pf) is the electron energy, 
counted from the Fermi energy, T, R ' A (e) is the electron 
retarded or advanced self energy, and v-p is the Fermi ve- 
locity. Performing the integration over momentum p, we 
obtain 



1 



de 



3 J ImS^(e)4Tcosh 2 £/2T 



(20) 



with v being the density of states of conduction electrons 
per one spin orientation. 

We assume that the electron self energy part contains 
two components: 



ImE A (e) 



1 

2t7 



3 ImT(e). 



(21) 



The first term in Eq. J3TJ, l/2r e , represents the effect 
of elastic scattering off non-magnetic impurities with r e 
being the mean free elastic time. The second term, 
n s ImT(e), represents the effect of scattering of electrons 
with energy e off magnetic impurities. Here the scatter- 
ing off a particular magnetic impurity is characterized by 
the T— matrix T(e); the self energy E^e) contains T(e), 
averaged over various impurities. 

Using the simple relation p = l/er between the conduc- 
tivity a and the resistivity p and Eqs. (|2*0|l and (|2~TT) . we 
represent the resistivity as a sum of two terms 



P 



(22) 



The first term is the resistivity of a metal without mag- 
netic impurities (n s = 0), which is produced by elastic 
scattering off non-magnetic impurities: 



1 



Pc 



er c = 2e 2 vD : 



(23) 



where D = VpT e /3 is the diffusion coefficient. The sec- 
ond term is the contribution to the resistivity due to the 
scattering off magnetic impurities: 



3n s 



ImT(e) 



de 



4Tcosh 2 e/2T 



(24) 



The scattering T— matrix in Eq. I|24(l has different struc- 
ture in the limits of high (T > 7k) and low (T < 7k) 
temperatures. We study these two limits below. 



A. High concentration of magnetic impurities, 

T sg > T K 

At high temperature T 7k the scattering of elec- 
trons off magnetic impurities is described by the Born 
approximation with the exchange constant renormalized 
according to Eq. (J2J. In this case the T— matrix is (see 
Appendix lB|l : 



ImT(e) = vvj 2 J K(u))j^ 



s £ / T dio 



+ e (e-u)/T 27r ' 



(25) 



where K(lj) is the Fourier transform of the spin corre- 
lation function, defined by Eq. Substituting ImT'(e) 
into Eq. fZQ. we obtain the Kondo contribution to the 
Drude conductivity: 



3n s J 2 
2e 2 v 2 



K(lu) 



dw 



Tl 



j/t- 



(26) 



We emphasize that Eq. I|26|l is valid if the distance be- 
tween impurities is much larger than the Fermi wave- 
length. At the same time, Eq. H26JI describes the resis- 
tivity in metals with arbitrary structure and strength of 
interaction between magnetic impurities^ We perform 
further calculations using the spin correlation function 
K(uj) given by Eq. Ijl5|l . which was derived within the 
virial expansion. 

Using K{oj), calculated within the virial expansion 
Eq. (|15f) . and the Kondo renormalized exchange constant 
J, see Eq. ©, we obtain 



127m s 5(5 + 1) 
e 2 vlv 2 \n 2 (T/T K ) 



1 - a s - 



Tsg 
T 



(27) 



Here numerical coefficient as is given by the following 
integral: 



+ 00/ 



as = 



1 - 



2S , 
,7.7' 



ej,)e~y^ (2J+1)Ajj, \ dy 



1 



-viej-ej,) S{S+\)Z{y) / y 2 



(28) 

where ej, Z(y) and Ajji are defined by Eqs. I|ll|) . I|13a|) 
and (|13b|) respectively. We emphasize that the integral 
in Eq. I|28fl converges near y — 0. The values of as are 
presented in Tabled 

Equation 127|l is similar to the results of Refs. I8HT7I 
Unlike Ref. \VA Eq. I|27() takes into account the Kondo 
renormalization of the exchange constant, Eq. @, and 
of the RKKY interaction, Eq. (|1L)|). The advantage of 
Eq. H27|) in comparison with Ref. |8| is in a consistent def- 
inition of the energy scale T sg and in an accurate proce- 
dure of the averaging over states of magnetic impurities, 
which allowed us to calculate the numerical factor as- 

In metals with low Kondo temperature, 7k ^ T S g, the 
competition between the Kondo effect and the effect of 
RKKY interaction results in a maximum of the resistivity 
as a function of temperature. If 7k ^ TU, then the 
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TABLE I: Values of the numerical coefficient as, Eq. I28L 
for several values of S. 



s 


1/2 


1 


3/2 


2 


5/2 


3 


7/2 


a s 


1.10 


1.52 


1.85 


2.11 


2.33 


2.52 


2.69 



maximum occurs at 



1 ~~2 sg 



T 
lnX 
T K 



(29) 



see also Ref. ItI IsIoL We notice that temperature T : is 
within the region of applicability of the virial expansion, 
used in the derivation of Eq. 1)27(1 . 

As temperature T approaches and crosses T sg , intrin- 
sic random magnetic field develops, and the renormal- 
ization of the exchange constant is stopped at J ~ 
2/(i/ln(T sg /Tk)). Simultaneously, the virial expansion 
breaks down, and the collective modes of spin system 
have to be considered in the derivation of the tempera- 
ture dependence of the resistivity. 

At present, there is no theo ry of m etallic spin glasses, 
which would provide us with K{ui) in a broad range of 
uj. Thus, the explicit form of Apk(T), see Eq. (|26H . is 
not known. We expect that A/?k continues to decrease 
with temperature T decreasing, as the dynamics of the 
local magnetic moments gets progressively suppressed at 
lower temperatures. In the mean field picture, each of 
the spins at T = is subject to a finite field, so that such 
dynamics is fully suppressed. In this case, the limiting 
value of A/ok at T = can be estimated as 



Ap K (T = 0) 



127T7J, 



s 2 



e 2 v^ 2 In (T sg /Tk) 



(30) 



(The quenching of the spin flips leads to the replacement 
of S(S + 1) factor, see Eq. (27}. by the factor S 2 here.) 
The picture leading to Eq. (|30|l can not be valid for the 
"tightest" pairs of magnetic moments with the charac- 
teristic interaction energy significantly exceeding T sg . It 
is not clear to us at the moment even in which direction 
the estimate (|27f) changes due to the deviations from the 
mean field description. 



B. Low concentration of magnetic impurities, 

T sg « T K 

At high temperature (T ^> 7k), the resistivity is de- 
scribed by Eq. J27J even in the limit 7k 3> Tsg . Neverthe- 
less, the effect of interaction between magnetic impurities 
is small due to the factor T sg /T in Eq. 1)2 7|l . and the max- 
imum of the resistivity does not occur at T > 7k- We 
show that at temperatures T < 7k the resistivity also 
monotonically increases as temperature decreases. 

For this purpose we use the following form of the imag- 
inary part of the scattering T— matrix: 23 



Here the first term is the contribution to the T— matrix 
in the unitary limit and T(e) is the T— matrix, written 
for the residual interaction between conduction electrons 
and magnetic impurities. The form of T(e) differs for spin 
impurity S — 1/2 (exactly screened magnetic impurities) 
and for S > 1/2 (underscreened magnetic impurities). 23 
We consider the two cases in more details below. 



1. Spin S =1/2 impurities 

In this case the majority of magnetic impurities are 
completely screened and the RKKY type interaction be- 
tween them is absent (only a small part oc T sg /Tk of mag- 
netic impurities form coupled states with binding energy 
exceeding 7k). Therefore the contribution to the resis- 
tivity is determined by the T— matrix of a single mag- 
netic impurity, which is well studied at T <C Tki^&Si 
We present the results for convenience. The scattering 
matrix for the residual interaction between conduction 
electrons and magnetic impurities has the form: 



Imf(e) 



1 9tt 
^87f 



{is 2 



TT 2 T 2 



(32) 



Substituting Eqs. <|HD and into Eq. and per- 
forming integration over energy e, we obtainiS 3 . 



Ap K = Api. 



1 ^Tt 2 



Here the factor 



A Pu = -- 



(33) 



(34) 



ImT(e) 



Imf(e). 



(31) 



corresponds to the unitary contribution to the resistivity 
at T — 0. According to Eq. (|33|l at finite temperature 
the resistivity contains corrections to Eq. H34J) ■ which are 
proportional to T 2 /T^. 

We also notice that the coefficients in Eq. I|33|l con- 
tain small corrections T sg /TK due to magnetic impurities, 
which form coupled states with binding energies exceed- 
ing the Kondo temperature. Two possibilities exist: i) if 
the coupled state is a singlet, these impurities do not af- 
fect electron transport; ii) if the coupled state is a triplet, 
it becomes screened and again leads to a T 2 dependence 
of the resistivity on temperature. 25 We emphasize that 
the ratio of the number of such impurities to the total 
number of magnetic impurities is small as Tsg/Tk. 

As the impurity concentration n s increases, and conse- 
quently T sg increases, the system of magnetic impurities 
with spin 5 = 1/2 undergoes a quantum phase transi- 
tion^ to a spin glass state. 



2. Impurities with S > 1/2 

As was shown in ref. |24j , at T <C 7k the residual cou- 
pling of conduction electrons with magnetic impurities is 
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described by the exchange Hamiltonian Eq. (JTJ with ef- 
fective impurity spin, S = S —1/2, and the renormalized 
exchange constant 



J = 



v\nT K /T' 



(35) 



This coupling results in the RKKY-like interaction be- 
tween magnetic impurities, which may be written in the 
form of the Hamiltonian given by Eq. (7J with the ef- 
fective spin S and the strength V(R) of the RKKY in- 
teraction defined as a solution of the following equation: 



The collective modes existing in the spin glass state 
would result, however, in deviations from Eq. H4U|) . 

To summarize Section ITTT1 the Kondo contribution to 
the Drude resistivity is a non-monotonic function of tem- 
perature only in the case of relatively high concentration 
of the magnetic impurities, i.e., at T sg > Tk- In the 
opposite case, this contribution increases monotonically 
with the decrease of temperature, at any value of the spin 
S of magnetic impurities. 



IV. PHASE RELAXATION RATE 



V(R) = 



ttvR 3 ln 2 [T K /U(I?)]' 



(36) 



Using Eq. I|36[) we estimate the typical value of the inter- 
action between magnetic impurities. Within logarithmic 
accuracy we have 



2n s 



irv In {T-^v /n s ) 



(37) 



cf. Eq. (TO). 

The imaginary part of the T(e)— matrix is given by 
Eq. (EH) with the impurity spin S = S — 1/2, exchange 
constant J given by Eq. I|35|) . and the typical value of 
the RKKY interaction between impurities T sg given by 
Eq. l|37)l. Substituting ImT(e) from Eq. with the 
modified parameters into Eq. l|3T|) and using Eqs. lt27)|) 
and lO, 

we obtain 



Ap K = Apt. 



1 - 4tt 



: S 2 ~l/4 
In 2 T K /T 



1 n Ts s 

1 - as-l/2y 



Unlike the previously analyzed case of high impurity 
concentration, T sg 3> Tk, here the temperature depen- 
dence of the resistivity remains monotonic. The lead- 
ing contribution to the (negative) derivative dp/dT in 
the temperature interval T* < T < Tk comes from the 
Kondo renormalization of the exchange constant; here 
the characteristic temperature T* is 



2 ss T sg 



(39) 



Below T* , the derivative dp/dT is determined by the in- 
teraction between magnetic impurities. 

As temperature approaches T sg , the virial expansion 
ceases to be valid, and the system may attain a spin glass 
state. Unlike the case of higher impurity concentration 
(T sg ^> Tk), here we expect the dependence p(T) to level 
off at T < T se .. This difference stems from the behavior of 
' i i scattering matrix, see Eq. I|31|) , in the vicinity of the 
unitary limit. If each of the local moments is quenched 
individually, then the saturation would occur at 



l-4r 



(S-l/2) 2 



hi Tx/T sg 



(40) 



Weak magnetic fields suppress the interference contri- 
bution to the conductivity, the difference 



Acr wl = a(B = 0) - o(B > B a ) 



(41) 



in the conductivity at zero magnetic field B — and at 
sufficiently strong magnetic field B S> B is called the 
weak localization correction to the conductivity. The 
characteristic value B of magnetic field, which sup- 
presses the weak localization is 



_ he / 1 
° ~ e V Dr^A' 



(42) 



where A is the cross-section area of the wire. Throughout 
this Section, we assume that the scattering off magnetic 
impurities dominates over all other mechanisms of the 
electron phase relaxation. 

We calculate Act w i of a metal with magnetic impurities 
and strong spin orbit coupling. In this case only the 
singlet component C s of the Cooperon remains finite^ all 
other Cooperon components are suppressed by the spin- 
orbit interaction. The weak (anti)localization correction 
to the conductivity is given byi 27 i 28 



Acr„ 



e 2 D f dede\de2 d d q 



c: 



:,£i 
-2,£ 



(2tt) 3 (27r) d 2Tcosh 2 e/2T 



(43) 



Neglecting the RKKY interaction between magnetic 
impurities, we have the following expression for the 
Cooperon: 



4tt 2 <5(£i - £ 2 )S(e[ - 4) 



C *°{e[J 2 ;q ) IV ! /I.-, - :',) f 2/-/ 



(44) 



where l/r s is the electron scattering rate off magnetic 
impurities and is defined by Eq. JHJ). 18 In calculating the 
WL correction, we assume that the rate 1/t s is higher 
than the Korringa relaxation rate of the magnetic im- 
purities, see Ref. for further discussion. Performing 
integration over the momentum q in a one dimensional 
case, d = 1 in Eq. (|4^|) . which adequately describes wires 
of the cross-sectional area A Dr., we obtain 



A<7. 



(0) 
wl 



e 

2irh 



2 



(45) 
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We study the effect of interaction between spins of 
magnetic impurities on the weak localization correction 
to the conductivity. As we have already mentioned, the 
interaction between magnetic impurities lifts the degen- 
eracy of impurity spin states, and therefore it is remi- 
niscent to the Zeeman effect of external magnetic field. 
The Zeeman splitting of spin states of magnetic impuri- 
ties affects the WL correction to the conductivityjiSiiii if 
the splitting is larger than either temperature T or phase 
relaxation rate l/r s . Similarly, the RKKY interaction be- 
tween two impurities starts to affect the WL correction 
if the RKKY interaction strength exceeds T or l/r s . We 
calculate the WL correction at temperatures T higher 
than T sg (or T sg ), so that only a small number of mag- 
netic impurity pairs satisfy this condition, and therefore 
the virial expansion is applicable. 

We notice that because of the RKKY interaction 
between impurity spins, the scattering processes may 
change electron energy and, particularly, may switch the 
position of the Cooperon poles in energy plane with re- 
spect to the real axis. These processes result in mixing 
of the Cooperon components C + and C~ which have dif- 
ferent analyticity, see e.g. Eq. (|44|l . The full equation for 
the Cooperon is 



• fc lJ £ 2 



• £ l! fc 3 



£3) £4 
-3; fc 4 



£4,^2 
£4, e 2 



(46) 



The diagonal elements of the 2x2 matrix C s are Cf . The 
matrix 



C s ,o(-,q) = 



C 



s.0\ 









Co (■,<*), 



(47) 



is the Cooperon to the zeroth order in the RKKY inter- 
action, see Eq. I|44|l . while the self energy 



" I Ef (■) Ef (•) 



(48) 



contains the higher-order RKKY contributions. 

To evaluate the first term of the virial expansion, we 
may account for £ s by the first-order iteration of the 
solution of Eq. f4*6"|l . The self energy £ must be calculated 
up to the first order in the RKKY interaction. In fact, it 
is sufficient to evaluate the upper diagonal element £"(•) 
of the matrix E s (-) and write the Cooperon as: 



1 



1 



X " " \e 2 , ej Dq 2 + 2/r s - i( Sl - e) 
with (see Appendix IP) 



Dq 2 + 2/r s Dq 2 + 2/t s - i{e - e 2 ) 

E,E 1 \ 1 



, (49) 



• 1 \ ,_o,... ? , / ^ (e-/ T + l)dK 2 (cj) 



e (e-u)/T + 1 

x {8{ei +lj — e)5(e — E2 — to) + 6(si - e)S(e ~ £2)} • 



Then, the weak localization correction to the conductiv- 
ity is obtained by substituting this expression for Cf (•; q) 
into Eq. (|43[1 and performing integration over momentum 
q and energies. As we discussed in the previous Section, 
details of the structure of the spin correlation function 
5K2{uS) in Eq. (|50|l depend on the relation between tem- 
perature T, Kondo temperature Xk and the typical en- 
ergy of interaction between impurities T sg . Some of these 
limits are discussed below. 



High concentration of magnetic impurities, 



At temperatures T ^> Tk the scattering of electrons 
off magnetic impurities is described by the Born ap- 
proximation with the renormalized exchange constant, 
Eq. J2Jl. We substitute SK 2 (uj) (the first order term in 
the RKKY interaction) from Eq. lfH>|l into Eq. (|5T)|l , and 
using Eqs. (|43|) and l]49fl. we obtain the weak localization 
correction to the conductivity [see Appendix iDl for more 
details]. We distinguish three temperature domains for 
the WL correction to the conductivity Aer w i. 

In the highest of the three domains, T ^S> 2/t s , the 
weak localization correction has the form 



Acr„ 



e 

2irh 



Dt s 
2 



tt(45'+1)(45' + 3) 
120(25+1) 



with the second term in the parentheses coming from the 
RKKY interaction. The use of Eq. for T sg and of 
the estimate for the Kondo-renormalized electron spin 
relaxation rate, 



1 



87m s S(S + 1) 
v In 2 T/T K 



allows us to estimate T sff r s as 



ln 2 (T/T K ) 



47T 2 5(5+l)l n 2 (z;Fn V3 /TK) 



(52) 



(53) 



(50) 



We see now that the correction due to the RKKY in- 
teraction only weakly depends on temperature, and is 
numerically small. 

It is curious to notice that the second term in Eq. (|5 1|) 
is almost independent of n s . This term takes into ac- 
count the fact that the contribution to the phase relax- 
ation rate is suppressed, if a scattering process results 
in energy exchange larger than l/r s . The reduction of 
the phase relaxation rate leads to the enhancement of 
the weak localization correction to the conductivity, as 
shown in Eq. I|51(l . The number of impurities with the 
splitting of energy states larger than 1 /r s constitute only 
T sg T s part of the total number of magnetic impurities. 
We emphasize that the accidental numerical smallness of 
the RKKY-induced correction oc T sg r s justifies the use of 
the conventional theory— of the weak localization in the 
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presence of magnetic impurities in the considered tem- 
perature domain. 

In the second temperature domain T sg < T < 1/t s , the 
WL correction equals 



Afiwi 



e 

2nh 



Dt s 
2 



1 



as_Tsg 
2 T 



In- 



T 
Ik 



1 



2 T 



(54) 



Here the numerical factor as is defined in Eq. I|28|l . and 
T sg is given by Eqs. ifTUjl . The result shown in Eq. lf5^|) 
has a similar structure to the expression for the resis- 
tivity correction Eq. Q27f) . The dependence Aer w i vs. T 
has a minimum at temperature T* defined in Eq. (|29H . 
This minimum results from the competition between two 
opposite trends: with the reduction of temperature, r s 
gets shorter, see Eq. (|52|) . while the stronger-bound im- 
purity pairs stop affecting Acr w i. Note that due to the 
relation between T sg and l/r s , see Eq. (|53|l . the second 
temperature domain is rather wide. 

The third temperature domain corresponds to the spin 
glass state of the magnetic impurities. With temperature 
decreasing to T sg , the virial correction becomes large and 
Eq. I|54|) is no longer applicable. At such temperature a 
spin glass transition is expected. Below the transition, 
Aer w i is still determined by the spin correlation function, 
Eq. ©. Similar to the discussion of the resistivity in 
Section IlIII we expect a monotonic increase and satura- 
tion of the WL correction. The limiting value of Act w i 
at T = was estimated in Ref. |2!j where quenching of 
the dynamics of each of the local moments was assumed. 
Deviations from such a simple picture of the spin glass 
state would result in a different value of Act w i(T = 0). 



B. Low concentration of magnetic impurities, 

T sg < T K 

At T ^ 7k, the weak localization correction to the 
conductivity is still given by Eq. I|54f) . However, now the 
effect of the RKKY interaction on Aer w i is small, and the 
WL correction to the conductivity decreases monotoni- 
cally with the decrease of temperature. 

At temperature T <C 7k the scattering off a single im- 
purity approaches the unitary limit. The potential scat- 
tering characterizing the unitary limit, does not destroy 
phase coherence and thus does not affect the WL cor- 
rection. Therefore, at T <C 7k only small deviations 
from the unitary limit determine Act w i. In this section 
we show that Act w i increases monotonically as tempera- 
ture decreases in the domain T <C 7k; the details of the 
temperature dependence are different for 5=1/2 and 
S > 1/2. 

Comparing the behavior of WL correction in the do- 
mains of low and high temperatures, we conclude that 
the correction must have a minimum at T ~ 7k, assum- 
ing that the scattering off magnetic impurities dominates 
the electron phase relaxation. 



1. Spin S= 1/2 

At T = 0, the spins of magnetic impurities are com- 
pletely screened and do not contribute to the phase 
relaxation of the conduction electrons. At finite but 
small temperatures, T <C 7k, the residual local electron- 
electron interaction facilitated by local moments leads to 
the electron relaxation which affects the Cooperon pole: 



,1/2 



£,£l_ 



4ir 2 S(e - ei)S(e - e 2 ) 
Dq 2 -i(ei -e) +f(e)' 



Here the relaxation rate 



T(s) 



9tt n s 3e 2 + n 2 T 2 
Y~ f 2 



(55) 



(56) 



Because the interaction responsible for the relaxation is 
local, the typical energy transferred in a scattering event 
is Ae - T, and therefore f (e)/Ae < 7 1 T sg /T K < I. Un- 
der these conditions, the Cooperon relaxation rate F(e) 
is just twice the one-electron relaxation rate. 

Substituting C + (-,q) from Eq. JHSJ into Eq. (15}. we 
obtain 



Act- 



wl 




de 



r(e) 4Tcosli e/2T 



(57) 



According to Eq. H57JI . in the absence of other phase 
relaxation mechanisms, the weak localization correction 
would vary as l/T at T < T K : 



Act- 



wl 



e 2 [vDT K 
0.022— W -. 

n V n s T 



(58) 



We notice, that similar to Eq. I|33|l . there are small 
corrections of the order of T sg /7k to the numerical co- 
efficient in Eq. (J5SJ. The corrections originate from the 
rare configurations of "tight" pairs of magnetic impu- 
rities, which form singlet or triplet states with binding 
energy exceeding the Kondo temperature. 



2. Spin S > 1/2 

The residual coupling between conduction electrons 
and magnetic impurities with S > 1/2 is still described 
by the exchange Hamiltonian Eq. with the reduced 
spin operator S = S— 1/2 and the renormalized exchange 
constant J , see Eq. (|35|l . In this case the electron scat- 
tering rate is 



1 



87m s S 2 - 1/4 
v In 2 T K /T' 



(59) 



The coupling J also results in the RKKY interaction 
between the partially screened local moments, which is 
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represented by the Hamiltonian Eq. J7J with S replaced 

by S. The strength of the RKKY interaction is deter- 
mined by the self-consistent equation (|36[1 . 

To calculate the weak localization correction to the 
conductivity, we use Eqs. I)43|) and H49|) with the self en- 
ergy E rr in the form of Eq. (|50|l . The spin correlation 



function SK2{uj) in Eq. l|BT]| describes correlations of S 
spins with the appropriately replaced exchange constant 
Eq. i|35|) and the typical value T sg of the RKKY potential, 

Eq. J37|). 

Similar to Section llll Al we can define three domains 
for the temperature dependence of the weak localization 
correction. At high temperature T > 2/f s , we obtain 



Act- 



wl 



27T?i V 2 



+ 2405 ss s 



(60) 



Equation l|6l7|) is a counterpart of Eq. (|ST|) . 

At lower temperatures, T sg < T < 1/t s , we obtain 
[compare to Eq. (|54"Jl ] : 



Aer w i = 



2tt^ 



Di 



&S-1/2 Tsg 

2 T~ 



(61) 



We notice that at T <C Tk as temperature decreases 
the weak localization correction to the conductivity in- 
creases. We conclude that both the Kondo effect and the 
effect of interaction between magnetic impurities reduce 
the phase relaxation rate as temperature decreases. At 
T* < T < Xk the temperature dependence of the weak 
localization correction to the conductivity is mainly de- 
termined by the Kondo effect, and at lower temperatures, 
T < T* it is determined by the interaction between mag- 
netic impurities. At temperature below T sg a spin glass 
state may appear. Similar to the behavior of Apx(T), 
the weak localization correction increases monotonically 
with the decrease of temperature, and should saturate at 
T ^ 0. 

To summarize Section IIVI the weak localization cor- 
rection to the conductivity is a non-monotonic function 
of temperature. At relatively high concentration of the 
magnetic impurities, i.e. at T sg > Tk, the positions of 
minimum in Act w i and maximum in Apx roughly coin- 
cide, see Eqs. lj2*7|) and In the opposite case, the 
minimum in Aer w i occurs at T ~ Tk- 



V. ENERGY RELAXATION RATE 

Free magnetic impurities are an intermediary for 
electron-electron scattering with small energy transfer A 
We show that the RKKY interaction between magnetic 
impurities leads to the electron energy relaxation as a 
result of a single electron scattering off a magnetic impu- 
rity. Indeed, if the impurity interacts with one or more 
of its neighbors, a scattering process is accompanied by 
the energy exchange between conduction electrons and 
magnetic impurities. 



In this Section we apply the virial expansion method to 
derive kinetic equations for the non-equilibrium distribu- 
tion function of electrons in a dilute magnetic alloy. The 
virial expansion is justified for processes with large com- 
pared to Xs g energy transfer from an electron to the sys- 
tem of localized moments; here T sg is the typical energy 
of interaction between magnetic impurities, see Eq. i|10fl . 
We assume that the spin-orbit interaction is strong, in 
which case the electron distribution is independent of 
spin orientation, / T (t,r,e fc ) = f L (t,r,e k ) = f{t,r,e k ). 

First, we consider electron scattering by magnetic im- 
purities belonging to a small-size pair. If the electron 
distribution function does not significantly vary on the 
length scale of the order of the pair size, then the corre- 
sponding scattering rate can be expressed in terms of the 
electron distribution function /(t, r, e k ) at the position r 
of the pair, 

r J k £,{t,r) = 2nJ 2 (2J +i)Ajj,5(e k - e k > + Ej - Ej,) 
x Pj(t, r, V)f(t,r,e k )(l - f(t,r, e k ,)). (62) 

Here Pj(t, r,V) is the distribution function for two mag- 
netic impurities over quantum states characterized by the 
total spin J of the pair, Ej — Vej, see Eq. lfTT|l . and Ajji 
is defined by Eq. (|13bjl . 

Having the rate T kk i (t), we can write the kinetic equa- 
tions for the distribution function of the pairs Pj(V). 
Performing summation over all initial (k,a) and final 
(&', a') states of a scattered electron as well as over the 
final states of the pair, we obtain the following equation: 



dPj(t,V) 
dt 



4^ 



2 J 



j' 



de k de k i 



( T JJ ' 

I 1 kk' 



1 k'k 



(we omit the position r of the pair in the argument of P). 
The normalization condition for Pj(t, V) has the form 



2S 
,7=0 



(64) 



In the stationary state, the distribution Pj(t,V) 
Pj(V) satisfies the equations 



P J+1 (V) = Pj(V) 



Jf(e)(l-f(e + Vej, J+1 ))de 
Jf(e + Ve JtJ+1 )(l-f(e))ds' 



(65) 



where we use the shorthand notation ej j> = ej — ej>, and 
ej is defined by Eq. (JTTJ. If the system of magnetic impu- 
rities and electrons are at equilibrium with temperature 
T, the solution of Eq. I|65|l is the Gibbs distribution: 



Pj{V) 



eM-Vej/T) 



E J (2J+l)exp(-y £j /T) 



(66) 



Next, we write the kinetic equation for the electron 
distribution function 



d 



d 2 



dt d 2 r 



f(t,r,e) = -I(t,e k ,V), 



(67) 
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where the electron collision integral has the form: 

J(t, s k , V)= v [ (Tj# - T&f ) dew (68) 
jj> J 

with Y^' given by Eq. 

We assume that the electron distribution function 
changes slowly with the coordinate r, so that the col- 
lision integral X may be averaged over a small volume of 
the metal, where /(t, r,£fc) does not change much, but 



which contains many magnetic impurities. In this case 
we can perform averaging of the collision integral over 
the RKKY potential according to: 

(T(t,e)> = J ^l(t,e,V). (69) 

Substituting Eq. (|62|l into Eq. (|68|l and performing aver- 
aging according to Eq. 1|69[) . we obtain 



, , NV 4T se ^ 2J+1 , , . f dE , 



I L:ll^[t,^)-(l-f e )f^ E Pj,(t,-^)\, (70) 



r 



where we use notations ejj> = ej — ej/ and f e = f(t, r, e). 
In the derivation of the collision integral, we tacitly as- 
sumed that the transferred energy E exceeds the width 
of spin states, given by the Korringa relaxation rate, 
h/T T ozTv 2 J 2 {T). 

The system of equations <|63[) and (|67fl , with the scat- 
tering rates given by Eq. Q62[l. the electron collision in- 
tegral replaced by its average, Eq. (fTPf . and with the 
initial conditions for f(t = 0, e, r), Pj(t — 0, V) and the 
boundary conditions for /(t, e, r), define completely the 
kinetics of electrons and spins. In the stationary case 
(df/dt = 0), one may use Eq. lffi5)l instead of Eqs. l|H2"|) 
and JHSJ- Note that the impurity average collision in- 
tegral, Eq. H70|). differs from the conventional collision 
integral for the electron-electron scattering. The l/E 2 
behavior of the kernel in Eq. (|70[l does not imply the 
scaling of the distribution function found in Ref. |2j 

The collision integral, Eq. (J70J), may be simplified for 
the electron energies e>T. As a result, we obtain the 
following kinetic equation for the distribution of "hot" 
(e ^> T) electrons: 



D 



d 2 f(e,x) 
dx 2 



= ±-Jx(E){f(e,x) 



dE 

f(s + E,x)}—. 

J-sg 

This equation is a version of a full kinetic equation, de- 
fined by Eqs. (|63|l - i|70|) . which may be used for analysis of 
the high energy tail of the electron distribution function 
in a metal with magnetic impurities. The kernel x(E) in 
Eq. (|71|l is asymmetric with respect to the energy trans- 
fer E and for |J5| 3> T sg has the following form 



X(E) = 



1ST 2 



„E/T( e E/T 



1) 



3 E 2 (3e B/T + l)(3 + e s/T )' 



(72) 



According to Eq. I|72() the probability for an electron to 
scatter with energy gain {E < 0) is exponentially small at 
l-El ^> T, while the rate of scattering with an energy loss 
E ^> T scales as a power of transferred energy, cx T 2 g /E 2 . 



We notice that because of the relatively slow decay of 
x(E) with energy E > 0, the relaxation of the number 
of non-equilibrium electrons occurs differently from the 
relaxation of their energy. We illustrate the energy trans- 
fer from electrons to the system of magnetic impurities 
by considering the following model problem. Assume, 
that initially the system of electrons and magnetic im- 
purities is in equilibrium at temperature T, and then 
instantaneously the electron subsystem is brought out of 
equilibrium, so that the new distribution function is char- 
acterized by small deviation Sf(e) from the equilibrium. 
The excess electron energy W per unit volume may be 
defined as W = 2u j e6f(e)de. The energy W will de- 
crease in time as the result of the energy redistribution 
between electron and impurity subsystems. Eventually, 
a new equilibrium with new temperature will establish. 
We calculate the reduction of the electron energy at the 
initial moment. The result is 



dW 
~dT 



16T f 



3t s 



dE 



Z 2E/T _ 1 



E (3e B/T + l)(3 
v8f{e) sinhe/T 
cosh^/T + coshe/T 



5/T) 



(73) 



de. 



If the distribution of electrons at the initial moment was 
peaked near energy e ^ T, e.g. v5f(e) = aS(e — e ), 
then the estimate of the energy reduction rate is dW/ dt = 
aT sg /T s ln£o/T. The characteristic collision rate is ~ 
1/t s and the typical energy transferred in a collision is 
T sg . However, the range of the transferred energies is 
broad enough to result in a logarithmic dependence on 
e /T. 

In conclusion we emphasize that Eq. (|72|l is derived 
within the virial expansion and is valid for large energy 
transfer. We expect that beyond the virial expansion the 
kernel x(E) remains to be a function of E/T sg and T/T sg : 



X(E) = T 



E T 

J- sg J-sg 



(74) 
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b 
< 



T* 



T 



FIG. 1: Schematic picture of the Kondo contribution to the 
resistivity ApK and of the weak localization correction Aa w i 
to the conductivity in the limit T sg S> Ik- Both ApK and 
Actwi are expected to saturate as T approaches T sg . 



Function T characterizes the excitation spectrum of a 
system of magnetic impurities at the energy scales rel- 
evant for the kinetics of conduction electrons. The 
study of its properties may provide important informa- 
tion about formation of spin glass states in metals with 
magnetic impurities. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we considered the effect of the RKKY in- 
teraction between Kondo impurities in a metal on kinetic 
properties of conduction electrons. Specifically, we evalu- 
ated the effect of interacting magnetic impurities on: the 
momentum relaxation rate and the corresponding contri- 
bution to the Drude resistivity of a metal Apx, the phase 
relaxation rate as denned by the weak localization cor- 
rection Aa w i to the Drude conductivity, see Eq. l)3T)l. and 
the energy relaxation rate, which determines the relax- 
ation of non-equilibrium electrons injected into a metal. 

The overall temperature dependence of the momentum 
and phase relaxation rates differs for the cases of strong 
and weak RKKY interaction between the magnetic im- 
purities. 

If the interaction T sg between the impurities sepa- 

— 1/3 

rated by a typical distance n s is strong, T sg 3> Tk, 
then the momentum and phase relaxation rates are non- 
monotonic functions of temperature T, with the maxima 
at T ~ T*, see Eq. JUJ [here T K is the Kondo temper- 
ature for a single magnetic impurity, n s is the concen- 
tration of these impurities, and energy T sg is defined in 
Eq. ijTUI) ]. Therefore, if T sg » Tk then one expects a 
maximum of the Kondo contribution to the resistivity 
and a minimum of the weak localization correction Acr w i 
at T ~ T* , see Fig. The positions of these extrema 
shift towards lower temperatures with the decreasing con- 
centration of the magnetic impurities. At lower temper- 
atures T < T sg when the spin glass state is formed both 



< 











i . 2S 2 -l/4 

- . OC 1 — 47T — tj -, 

-N In T/Tk 


A<T W , , 











T 



T* 



T 



FIG. 2: Schematic picture of the Kondo contribution to the 
resistivity ApK and of the weak localization correction to the 
conductivity Aa w i in samples with low concentration (T sg <C 
Tk) of magnetic impurities with spin S > 1/2. The solid lines 
represent ApK and Aa w i for interacting impurities. Dashed 
lines show ApK and A<r w i as if the interaction were absent. 



the resistivity and the WL correction to the conductivity 
saturate. 

In the opposite case, T sg <C Tk, the momentum relax- 
ation rate increases monotonically with the decreasing T, 
and eventually saturates at T — 0, see the lower curve 
in Fig. |21 The saturation level depends on the value of 
impurity spin S. Thus the Kondo contribution to the 
Drude conductivity is a monotonic function of tempera- 
ture. The spin-induced contribution to the phase relax- 
ation rate, on the contrary, has a maximum at T ~ Tk- 
If this contribution dominates over all other mechanisms 
of the phase relaxation, then the weak localization cor- 
rection to the conductivity, Act w i, has a minimum at 
T ~ Tk, see the upper curve in Fig. [21 The details of the 
low-temperature increase of Act w i with the further reduc- 
tion of temperature in the region T < Tk depend on the 
level of Kondo screening. In the case of full screening 
(S = 1/2), the phase relaxation rate vanishes at T — > 0, 
and Acr w i(T) diverges, see Eq. (|58[l. If the screening is 
not complete (S > 1/2), then both the phase relaxation 
rate 1 /f s and the weak localization correction Acr w i sat- 
urate at some finite level, see Fig. EI Within the simplest 
model of the spin-glass state accepted in Ref. |2^, we find 



Act 



wl 



e 

2^h 



2 



1 



8tto s (S - 1/2) 5 
v In 2 T K /f sg 



(75) 



where the temperature T sg is defined in Eq. (|37J) . It is 
worth noting that the saturation occurs at temperature 
~ T sg , well below the Kondo temperature Tk- 

The considered limits of T sg > T K and T sg < T K , 
and the conjecture of Herta 2 ^ allow us to understand the 
evolution of the temperature dependence of ApK (T) and 
Act w i with the concentration of impurities n s . The po- 
sition of the maximum in Apx{T) shifts continuously 
towards T = with the decrease of n s ; it reaches T = 
at some finite value of n s , which corresponds to T s<r ~ Tk- 
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Note that such behavior occurs irrespective to the value 
of S. Formation of the spin glass with T sg <C Tk a t 
S > 1/2 does not result in the finite-temperature maxi- 
mum of the function Apx ■ The position of the minimum 
in A<7 w i also moves to lower values of T with the de- 
crease of n s . This shift, however, stops at T ~ 2k; thus 
the minimum occurs at a finite temperature even in the 
limit T sg <C 2k- (Once again, here we assume that t s is 
the shortest of the phase relaxation times.) 

Electron scattering off interacting magnetic impurities 
leads to the energy transfer from electron to the system 
of localized spins. The rate of collisions with a relatively 
large energy transfer E can be calculated by means of 
the virial expansion. The corresponding full system of 
kinetic equations for the electrons and spins is derived in 
Section[V] The collision rate with energy loss E at \E\ 3> 
max{T sg , 2k, T} scales with E as 6(E)T se /(t s E 2 ). This 
asymptote of the rate is not sensitive to the formation 
of the spin glass state. However the spin glass transition 
affects the electron energy relaxation for smaller energy 
transfers, \E\ < T sg . 

In most part, the data of existing experimental 
can be understood within the pre- 
sented here theoretical framework. The evolution of the 
temperature dependence of the resistivity with the con- 
centration of magnetic impurities was studied in dilute 
AuFe alloys. The investigated range of the magnetic 
impurit y (Fe) c oncentration n s covered by the data of 
Refs. lo MilUll lM3( j| .3l l is extremely broad, ranging from 
3.3 ppm, see Ref. [l(J, up to a few percent^ .. For concen- 
trations n s > 100 ppm, the function Apx(T) has a clear 
maximum^S. Its position T max moves to lower temper- 
atures as the impurity concentration n s decreases; the 
measured in Ref. dependence of T max on n s is super- 
linear, in a qualitative agreement with Eq. I|29|l . The 
extrapolation of the data of Ref. |f| to T max = yields 
the critical value of Fe c oncentr ation « 50 ppm for the 
AuFe alloy, see also Refs. ll0l3lL Finally, the observed^ - 
magnetic hysteresis of the resistivity at n s = 7 ppm may 
indicate formation of a spin glass even for such low im- 
purity concentrations, which is possible at S > 1/2, see 

Section inrpHi 

The weak localization correction to the conductivity in 
AuFe wires with low impurity concentration (n s = 7 ~ 
60 ppm) was studied in Refs. limn There is a proper 
correspondence between the dataii for the sample with 
n s = 60 ppm, and the data 10 for n s — 10.9 and 7.1 ppm. 
The values of l/r s found from the weak- localization mag- 
netoresistance, see Eqs. I|52|l. <|5^)l and J7SJ), scale roughly 
linearly with n s . The temperature dependence of l/r s for 
the investigated samples is also in accord with the the- 
ory. Namely, the phase relaxation rate exhibits a broad 
plateau at temperatures around 2k ~ 0.3 K (with the 
plateau valueii of l/r s w 6 x 10 10 s _1 for n s = 60 ppm). 
The plateau is followed by a decrease of this rate at lower 
temperatures. The saturation^ of l/r s (T) at T < 2k is 
compatible with the value of impurity spins S > 1/2. 

We also notice that the data of Ref. ^] for an AuFe 



alloy with n s = 15 ppm are in sharp contrast with 
other experimental data^i2iii*2i and with the expecta- 
tions supported by the presented theory. Indeed, Ref- 
erence [ll| reports the position of the Drude resistivity 
maximum at T max ss 30 mK for the n s = 15 ppm sam- 
ple, which is indistinguishable from the value of T max for 
the n s = 60 ppm sample in the same work. This stability 
of T max contradicts the dependence of T sg (n s ), expected 
from other experimental works^i and from theory, see 
Sec. 11111 Also, the four-fold difference of n s between the 
two samplesii resulted in a 100-fold decrease of the elec- 
tron phase relaxation rate l/r s . This drastic change of 
1/t s with n s contradicts both the quoted above measure- 
mentsi - and the theoretical estimates, see Eqs. (|52|l and 

Measurements of the energy relaxation in nano-wires 
of Au, Cu, and Ag revealed the effect of individual mag- 
netic impurities 2 , but there was no systematic study of 
the effect of RKKY interaction on the electron energy re- 
laxation. Measurements of the relaxation rates at energy 
transfers \E\ < T sg may provide information about the 
excitations in a spin glass, but we are not aware about 
such measurements as of yet. 
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APPENDIX A: MATRIX ELEMENTS Ajj, 

In this Appendix we present the expression for the 
Ajj* factors in Eq. I|12|) . First we notice that the con- 
venient form to calculate ( JM\S a \ J'M') is to use the 
basis of spin states for two independent spins: \m1m2), 
where m\ and 1112 are the spin components along some 
direction. We have 

r) a (JMJ'M') = (JM\S a \J'M') 

= Cmim 2 ( T7l l|'5'ct| TO l)Cm' 1 7Ti2' 

mi, 7712,^1 

where the Clebsh-Gordon coefficients are expressed in 
terms of the Wigner 3j— symbols as 

c- 2 = V27Tl(_ J M I I). (Al) 

The matrix element Ajji may be represented in terms 
of t] a : 

A "'= E 52vl(JMJ'M'). (A2) 

a—x,y,z M' 
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Reo; 



2mTm 



FIG. 3: The contour for calculation of the T— matrix, 
Eq. HB4I . in the Matsubara representation. The contributions 
from the dotted parts of the contour vanish. 



We find that 

(J + 1)(2S- J)(2S + J + 2) 
= jtt^ | ^ , (Ada) 



.4 



j, j -i 



4(2J+ 1) 
J(J+1)(2J+1) 

4(2J+1) 
J(2S- J + l)(25 + J+l) 



(A3b) 
(A3c) 



4(2J+1) 

and all other elements vanish. From Eqs. (|A3|) we verify 
explicitly that 

2S 

(2J + l)Ajj, = S(S + 1)(2S + l) 2 . (A4) 

J,J'=0 



APPENDIX B: ELECTRON T MATRIX. 

For completeness of the presentation, we show how the 
imaginary part of the electron T— matrix may be related 
to the spin correlation function, defined by Eq. J3J). In 
the Matsubara representation the T— matrix is given by 
the following expression: 

T(i£ n ) = J 2 TJ2 J pJ}3£(»Wm)0(*(E„ ~ W m ),p), 

(Bl) 

where K,{iu> m ) is the Fourier component at uj m = l-nrnT 
of the Matsubara spin correlator 



ic(T) = j2p^ {E( ~ Ei ' )T \(z\ s \z 



'\l 2 



(B2) 



[compare to Eqs. 10 and |0J for the real time spin cor- 
relator K(t)], and 



p) 



1 



(B3) 



is the electron Green's function at Matsubara frequency 
e n = n(2n + 1)T [compare to Eq. (|19|l ]. The integration 
over momentum in Eq. (|B 1 1) gives 



T(ie n ) = -iixvj' 2 £(ziL> m )sign(e n 



(B4) 



Next we perform the standard procedure of the analyt- 
ical continuation in Eq. (|B2|) . We replace the sum over 
discreet u> m by the integral over complex uj: 



f duj uj 
Tj2F(iw m )=i — coth — T{ 

J c ~ ■ ■ 



(B5) 



This procedure is valid for an arbitrary function an- 
alytic inside the contour C of integration. We choose C 
to be a circle of infinite radius with two cuts at Ivauj = 
and Imu = e n . Inside the contour of the integration, the 
function JC(uj) is analytic and only the poles of cotho;/2T 
contribute to the integral. Neglecting the contribution 
from the pieces of the circle, and keeping the contribu- 
tion along the cuts, we obtain 



T(iSn) 
TTvJ 2 



+ OG 

/ tanh— /C(ie„ + wj — 



+ 00 



(B6) 



+ signe„ / coth 



uj K{ui + iO) — fC(u> — iO) duj 



2T 



2tt 



The imaginary part of the T— matrix is given by 

T{e + iO)-T(e-iO) 
ImT(e) = i . (B7) 



We substitute Eq. (|B6|) in Eq. (|B7|I . take into account 
relation 



K{u>) = i 



K,{u) + iO) - K{u) - iO) 
e"/ T - 1 



(B8) 



between the real time spin correlator, Eq. J2J, and the 
Matsubara spin correlator, Eq. (|B2|) . and obtain Eq. 125|) . 



APPENDIX C: COOPERON SELF ENERGY 

The Cooperon self energy E rr (-) may be obtained as 

( zs^ \ 

a result of the analytical continuation of £ . ' . , I , 
written in the Matsubara representation: 



V:T : £i, e[\ _ (ex+iO, e[ + iO 
^ 1 e 2 , e'J [e 2 -iO, e' 2 -iO 



(CI) 



The Cooperon self energy is related to the scattering ma- 
trix S off a magnetic impurity for two electrons in a sin- 
glet state: 



S s (-) = 2-Kvn s S(-). 



(C2) 
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ie' 2 ,g' 2 ie 2^ 2 T ie' 2 ,g' 2 fe 2 .fc / ie' 2 , Q 2 

FIG. 4: The left two diagrams represent the self-energy cor- 
rections to the Cooperon and contain the single electron 
T— matrix. The diagram on the right is the vertex correction 
and is related to the two electron counterpart of a T— matrix. 



In Eq. i|C2l) S is averaged over various magnetic impu- 
rities. The singlet component iS s is related to the full 
matrix S of two electron scattering: 



'liirTm 



Imw ■ 



Rew 



(C3) 



Matrix <S, in its turn, contains the contribution from 
three diagrams, shown in Fig. 01 and has the following 
form: 



gSis'i ( l£ li t£ l\ _ y«i«i ( i£ lj %£ '\ 



^2 Vl£2, l £ 2 l ?2?2 \ t£ 2, i £ 2 

T(iei)-T{ie 2 ) af a a a a 

(C4) 

Here the second term contains the single electron 
T— matrix considered in Appendix [BJ Below we focus 
on the irreducible component of the scattering matrix of 
two electrons in a singlet state, represented by 



FIG. 5: The contour for calculation of the vertex correction 
V s to the Cooperon. The contributions from the dotted parts 
of the contour vanish. 



Next we perform the analytical continuation of V s . We 
consider the case e\ > and £2 < 0, see Eq. (|C1I) . Us- 
ing Eq. (|B5|> . we replace the sum over the Matsubara 
frequencies, u> rn , by the integral over a contour in the 
complex plane. In the present case the contour of the 
integration is shown in Fig. \5\ and contains three cuts: 
Imcj = £2, 0, e[. We notice, that the contour parts 
above the upper cut and below the lower cut do not con- 
tribute to the integral. The remaining parts of the con- 
tour along the cuts after the continuation ie[ — ► e\ + iO 
and ie 2 — > £2 — iO give 



E ( /T 



(C5) 



where ' (iuJ m ) is the scattering amplitude of electrons 
with initial spin state to the spin state <;, accompa- 
nied with the change of the state of magnetic impurities 
from £' to £. The summation over £ and £' runs over all 
possible states of the impurity spin and Z — e~ E z/ T . 
For scattering amplitudes in the Born approximation we 
have 



V,. 



+ OC 



x 2 



coth — [5K, 2 (w + iO) - SJC 2 (w - iO)] 



£ 2 - iO, el iO ) =J dLdS ( £l + u) ~ £ ^ £2 ^"-^ 

— OO 

- tanh ^ [8K2 (w + e[ + iO) - 6JC 2 (u + e 2 - iO)] } . 

(C8) 

This expression may be further simplified in the case 
when e[ = £ 2 with the help of Eq. (|B8f) . 



/«' = J(k\S\k'}&, 



(C6) 



We substitute Eq. ifHr)} into Eq. (fHsfl . and use Eq. (JHU 

to write down the singlet component of the vertex part 
of the Cooperon S matrix: 



y [ leu l£ i 
ie 2 , ie' 2 



1 7 2 0(-£ 1 £ 2 )0(-£' 1 £ 2 )r^ 5K 2 {iu m ) 



x5 



(C7) 



APPENDIX D: CALCULATIONS OF THE 
VIRIAL CORRECTION TO Acr wl 

To calculate the WL correction to the conductivity, 
we substitute Eq. I|49|) with the self energy defined by 
Eq. lj5U|) into Eq. (|4^|l . As a result we obtain 



Act 



wl 



Act 



(o) 

wl 



Act 



(l,o) 
wl 



Act 



(1,6) 
wl 



(Dl) 
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Here Aa^ is given by Eq. I|45|). and the second term is 



Act, 



(l,a) 
wl 



Act^ f 5K 2 (oj) uj 1 dw 



4 7 S(5+l)n-e-^ 7T 



Acr wi a s T sg 
2 T ' 



(D2a) 



with as defined by Eq. (JSSJ, see also Table |I] The third 
term in Eq. (|D1|) is 



As wr s — > function </?(w) vanishes. Therefore, at low 
temperature Tr s <C 1 only the second term Act,^' in 
Eq. IjDip remains. In the opposite limit, Tt s 3> 1, we use 
the following property of f(uj) 



tt-CILU = T, 

uj 2 32 



(1,6) _ Aa$ f 5K 2 {uj) uj ip{u) duj 



Aa 



wl 



2 J S(S + l)Tl~e-"/ T tt 



(D2b) 



and obtain 



where 



1 Im-v/1 + iu!T s /2 

U}(U)) = L ^^^^^^=. 

4 ^t sV /1 + ^ 2 t s 2 /4 



Act 



(1,6) 
wl 



tt (4S + l)(4S + 3) 



120 



25+1 



Act!$T ss t s , Tt s » 1. 

(D3) 
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